Mueller matrix polarimetry of plasmon resonant silver 
nano-rods: biomedical prospects 

Sayantan Ghosh*, Jalpa Soni*, Sudipta K. Bera, Ayan Banerjee and Nirmalya Ghosh 

Dept. of Physical Sciences, Indian Institute of Science Education and Research Kolkata, 
Mohanpur Campus, Mohanpur 741 252, Dist. Nadia, West Bengal, India 

ABSTRACT 

Fundamental understanding of the light-matter interaction in the context of nano-particles is immensely bene- 
fited by the study of geometry dependent tunable Localized Surface Plasmon Resonance (LSPR) and has been 
demonstrated to have potential applications in various areas of science. The polarization characteristics of LSPR 
in addition to spectroscopic tuning can be suitably exploited in such systems as contrast enhancement mech- 
anisms and control parameters. Such polarization characteristics like diattenuation and retardance have been 
studied here using a novel combination of MuUer-matrix polarimetry with the T-matrix matrix approach for 
silver nano-rods to show unprecedented control and sensitivity to local refractive index variations. The study 
carried out over various aspect ratios for a constant equal volume sphere radius shows the presence of longitu- 
dinal (dipolar and quadrupolar) and transverse (dipolar) resonances; arising due to differential contribution of 
polarizabilities in two directions. The overlap regions of these resonances and the resonances themselves exhibit 
enhanced retardance and diattenuation respectively. The spectral and amplitude tunability of these polarimetric 
parameters through the aspect ratios to span from the minimum to maximum ([0, 1] in the case of diattenuation 
and [0, tt] in the case of retardance) presents a novel result that could be used to tailor systems for study of 
biological media. On the other hand, the high sensitivity of diattenuation dip (caused by equal contribution of 
polarizabilities) could be possibly used for medium characterization and bio-sensing or bio-imaging studies. 
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1. INTRODUCTION 

Localized Surface Plasmon Resonances (LSPR) due to their geometry dependent tunability have been explored in 
the recent years in various contexts of biomedical applications like sensing, development of novel nano devices,^ 
probing bio-molecular interactions and their medical diagnostic and therapeutic applications^ just to name a 
few. Localized surface plasmon resonances while exhibiting potential in wide ranging applications have also 
been probed for understanding fundamental physical questions related to the interaction of light with nano- 
particles.'^'^ The localized field enhancement exhibited when resonance between the incident electromagnetic 
radiation and polarizability of the nano-structure is achieved; have been understood to stem from the intrinsic 
dependence of metal polarizability on the shape and size of the particle. The spectral position of such plasmonic 
resonances are also well known to have an inherent dependence on the ambient medium. Such properties of the 
LSPR have been exploited (as mentioned before) to tailor them for specific applications. 

The LSPR though traditionally investigated spectroscopically have recently been observed to show novel 
polarization properties^ which can be used as additional contrast mechanisms in the study of light matter in- 
teraction. The effectiveness of such polarization characteristics (for example) in bio-sensing applications can be 
understood to arise from the fact that the elastically scattered spectra from dielectric components in complex 
biological media generally can swamp the plasmonic spectral signatures. However, it must be noted that polar- 
ization parameters like diattenuation and retardance are spectrally independent for dielectric structures while 
enhanced for plasmonic structures. This realization makes the study of polarization characteristics of plasmon 
resonance of considerable interest in the context of possible applications in various areas of sensing. 
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Plasmon resonant silver nano-rods have been investigated in the present context of biomedical applications 
recently owing to their interesting LSPR properties, where, polarization characteristics like diattenuation (d, 
defined as the differential attenuation between two orthogonal polarizations), retardance {6, the phase difference 
between two orthogonal polarizations) and depolarization (A, the randomization of the polarized light after 
passing through a medium); of silver nano-rods are probed for different aspect ratios, s = d/l (where d and / 
are the diameter and length of the cylinder respectively) and ambient dielectric media. The T-matrix^ method 
for calculation of the Mueller matrices and a decomposition algorithm based on Lu and Chipman's method have 
been used as analytical and computational tools in this study. 

This paper is organized as follows: section 2 describes briefly the theoretical aspects of the work with a brief 
introduction to localized surface plasmons and the various analytical and computational tools used to study them. 
The extraction of the scattering matrices from the T-matrix method and their subsequent Mueller polarimetric 
study through a Lu and Chipman based approach is also described. Section 3 discusses the observations of 
the spectral characteristics of the LSPR for silver nano-rods along with the spectral polarimetric parameters 
diattenuation, retardance and depolarization as a function of their geometrical configurations. The dependence 
of the polarimetric parameters on the ambient medium refractive index is also investigated and a quantification 
of the sensitivity of diattenuation to the medium refractive index is provided. Finally a conclusion and possible 
future directions is discussed in section 4. 



2. THEORETICAL ASPECTS 



Surface plasmons can typically be described as collective oscillations of polarization charges generated at metal- 
dielectric interfaces, due to resonant illumination by electromagnetic waves and consequent conduction electron 
transport at the interface."^' The surface plasmons have been categorized as propagating (along the metal- 
dielectric interface) or localized (in sub- wavelength metal nano-particles). These resonances, owing to the po- 
larization response of the metal nano-particles are highly sensitive to the geometry, shape, size and ambient 
medium. The shape and size dependence of the resonances can be easily understood if one considers that the 
polarizability of a sub-wavelength metal nano-particle is a function of its volume and consequently shape and 
aspect ratio (e), which is the ratio of the lengths of the short axis and the long axis (e.g. for spheres, £ = 1). 

Several methods have been developed for the analytical (where tractable) and numerical calculations of 
electromagnetic scattering from small particles. It is worthwhile to note here that the scattering problem is 
analytically quite easy to deal with in case of rotationally symmetric particles like spheres, spheroids and ellip- 
soids. Notable among such methods are the Quasi-static or Electrostatic Approximation (ESA),^' Mie theory,^^ 
T-matrix method,^^'^^ Finite-Difference Time-Domain (FDTD) and Finite Element Method (FEM).!"^ 

In the electrostatic approximation, the sub- wavelength (a << A) size of the particle allows the approximation 
that the particle lies in a static field and hence the polarizability can be easily calculated. For metal nano-spheres, 
the dipolar polarizability a{Lu) and the scattering cross-section Cgcat can be calculated as^^ 
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Here, x — 27ra/A is the size parameter. It can easily be seen that for £i(w) — — 2£„(cj) is singularity point 
where, the Cscat will sow significant enhancement. This is the dipolar resonance condition for the dipolar sub- 
wavelength metal nano-sphere. Similarly, the condition for quadrupolar resonance is given £1(0-1) — — (3/2)£m(w). 
In general, for spheres, spheroids and ellipsoids, the dipolar polarizability given by 
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where, for spheres, a = b = c which implies that Li = L2 = = 1/3, for spheroids, a = b =/= c and for ellipsoids, 
a b ^ c. 

Although sub-wavelength particles can be easily studied using the quasi-static approximation, this method 
fails when subjected to particles whose sizes are comparable to that of the wavelength (r < A). Gustav Mie in 
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his pioneering work proposed the expansion of the incident and scattered electric fields in vector-spherical wave 
functions which could be used to solve the Maxwell's equations along with the boundary conditions (continuity of 
the electric and magnetic fields at the particle medium interface) exactly albeit only for spherical particles. ^"^^ 
The solution of the Maxwell's equations in this case take the form of spherical Bessel and Hankel functions. This 
method allows us to write the relation between the incident and scattered electric field components in the Jones 
formalism by Es = 3Ei where, the i^y and E± are parallel and perpendicular to the scattering plane respectively 
and J is the 2x2 amplitude scattering matrix or the Jones matrix which has a diagonal form for spheres. The 
elements of the amplitude scattering matrix J are represented by the scattering co-efficients a„ nd obtained 
from the Mie theory calculations. ^'^ The index n in the scattering coefficients represents the multipolar modes 
(e.g. n — 1,2 correspond to dipolar and quadrupolar modes respectively) and for each n, the corresponding a„ 
and bn can be interpreted as the transverse magnetic or TM (no radial magnetic field) and transverse electric or 
TE (no radial electric field) modes respectively.^" Since the Mie theory is applicable only for spheres, various 
extensions of this method to non-spherical particles like ellipsoids and cylinders have been proposed. However, 
the most prominent amongst them is the T-matrix method which calculates the scattering matrix for arbitrarily 
shaped particles by solving the Maxwell's equations while expanding the incident and scattered electric fields 
in terms of the vector spherical wave harmonics^^' and is computationally best suited to the calculation for 
rotationally symmetric particles.^' This T-matrix method analogous to the Mie theory yields the relation 
between the incident and the scattered electric fields's Stokes vectors {Si and Ss respectively) through a 4 x 4 
matrix called the scattering matrix or the Mueller matrix (named after Mueller for the Stokes-Mueller algebra), 
such that, Ss = MS^.^'^° The Mueller matrix by the virtue of its formulation contains all the polarization 
information about the scattering medium and the particle^® and can be explored using various Stokes-Mueller 
algebra based methods to extract and quantify polarization properties of scattering. ^^"^^ The polarization 
parameters thus obtained have been very important both in furthering the basic understanding of light-matter 
interaction^*^' and their applications to various front-line research areas like astronomy, meteorology,'^^ and 
biomedical research. ^^"^^ 

The study of the polarization characteristics of light matter interaction requires careful consideration of 
techniques that would provide the maximum information about the systems under study. For example, the 
Jones polarimetry involves coherent addition of amplitudes and phases of the electric field and is experimentally 
difficult to carry out. This technique also fails to account for the effects of partial polarization aspects. The 
Stokes-Mueller based technique is an intensity based technique and has emerged to be one of the most extensive 
polarimetric tools available to the community. As mentioned earlier, the incident and the scattered electric fields 
can be related by a 4 x 4 transfer matrix or Mueller matrix M which carries polarization information in both linear 
and circular bases. The elements of Mueller matrix M carry information like linear and circular depolarization 
(Al&Ac), retardance {5l^5c) and diattenuation {dL^dc) respectively. Depolarization is defined as the decrease 

. /Q2_{_f^2_|_y2 \ 

in the net degree of polarization (which is defined as DOP = — j ) due to the sample. Depolarization 

is given as A = 1 — DOP, where the incident EM-field is assumed to have a net DOP = 1. Depolarization can be 
understood to originate from multiple scattering processes in the sample. Retardance 5 is defined as the phase 
shift between two orthogonal polarization states. This property arises from the existence of polarization axes 
in axisymmetric particles like ellipsoids and cylinders. This phase difference for linear orthogonal polarization 
states like {H, V) and (P, M) is called the linear retardance 5l and for circular orthogonal polarization states 
like (L, R) is called as circular retardance 5c- Similarly, diattenuation d is defined as the differential attenuation 
of orthogonal polarization states while interacting with a sample and has the same origin as retardance. The 
differential attenuation of orthogonal linear (like {H, V) and (P, M)) and circular (like (L, R)) polarization states 
give linear (d^) and circular {dc) diattenuation respectively. 

Quantification and extraction of "lumped" polarization properties requires a mathematical decomposition 
of the Mueller matrix M, which can be done in various ways. One such method is the Polar Decomposition 
Method (PDM) following the prescription of Lu and Chipman.'^'^ The 4x4 Mueller matrix M can be decomposed 
into three "basis" matrices based on its polarization properties. This polar decomposition method sequentially 
extracts diattenuating M^, retarding Ts/lji and depolarization Ma matrices as 

M <^ Ma • • Mo, (3) 
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where represents the decomposition where, the decomposed basis matrices are defined as 




Pa ms J ' \0 mn 

Here, d — ■jj^[Mi2 M13 Mi4]-^ is the diattenuation vector, d — d/\d\ is the unit vector and mo = \/l — d^l + 
(1 — \/l — d'^)dd^ . The magnitude of the diattenuation vector \d\ gives the total diattenuation of the system 

and is expressed as d = V^i2 + ^13 + ^m- Similarly, Ps = ^^E^, where, P = jj^[M2i M31 M^i]'^ is the 
polarizance vector, ms and are 3x3 sub-matrices of Ma and Mj^ respectively and can be found by solving 
the eigenvalues of m' = ms ■ m^. The net depolarization can be calculated as A = 1 - M^^^^lnii^o < A < 1 
while the net retardance (i?) is given by i? = cos^^ ^ tr(Mfl) _ ^ r^j^^ linear retardance can be calculated from 

the Mfl as (5 = cos-i [^{Mr^22 + Mr^s^)^ + (M^^gs + 1/^,23)" - 1 • 

3. RESULTS AND DISCUSSION 

Silver nano-particles and nano-rods in particular have attracted the attention of the scientific community in 
the recent years in context to their varied applications. In this work, we consider silver nano-rods (circular 
cylinders) of equal-volume-sphere-radius (EVSR) 20nm and varying aspect ratios. Unless mentioned otherwise, 
the EVSR will be considered to be 20nm and the propagation direction of the incident electromagnetic field will 
be positive Z direction throughout the work. The refractive index data of Winsemius et al.^* has been used for 
the calculations. As has been mentioned before, the plasmonic signature of the metal nano-particle is affected 
by its orientation and the ratio of its short axis 2r [r being the radius) to its long axis I {I being the length 
of the nano-rod). Hence, we investigate the electromagnetic scattering and the polarization response of these 
silver nano-rods as function of their aspect ratio e = 2r/l. The detection plane is again the X ~ Z plane with 
the detection angle 6 = 0° signifying "on-axis" detection, i.e. the detection point lies on the propagation axis. 
6 = 90° detects the scattering at an angle perpendicular to the propagation direction; this can be either on the 
±X axes. In this study, we have chosen the scattering angle 6 = 45°. The scattering cross section Cgcat and the 
Mueller matrices have been obtained using the T-matrix codes written in Fortran?? by Mishchenko et al.^'^^'^^ 

The silver nano-rods by the virtue of a cylindrical geometry can be understood to have two polarization 
axes along the long and short axes of cylinder. For example, when the particle's long axis is aligned along the 
propagation direction of the incident radiation, a transverse axis is formed in the X — Y plane and a longitudinal 
axis is formed along the long axis of the particle. The particle size also plays a crucial role in the excitation of the 
plasmon modes. It is well known that due to the nature of the permittivity response of silver to electromagnetic 
radiation, the excitation of different orders of plasmon modes is limited. In fact, for spheres the minimum limit 
for quadrupolar excitation is 50nm.'* In rods, however, at an equal-volume-sphere-radius of 20nm, the aspect 
ratio 0.65 corresponds to a rod length of approximately 44.9nm and a rod diameter of 30nm, it is possible to 
excite quadrupolar resonances, albeit only in the transverse mode. The absence of the longitudinal quadrupolar 
mode can be understood easily from the consideration that the alternating positive-negative charge separations 
along the length of the cylinder cannot be sustained over such lengths. However, the longitudinal dipolar mode 
can be sustained over the length of the cylinder. 

Figure 1 shows the orientation averaged spectral variation of the extinction cross-section for nano-rods of 
EVSR 20nm at different aspect ratios. As expected, the quadrupolar plasmon mode is excited at wavelengths 
lower than the dipolar modes. The quadrupolar mode (excited at ^ 3?5 — 400nm) is transverse in nature and thus 
remains unaffected by the change in the aspect ratio which entails changes in the length of the nano-rods. The 
two dipolar modes, one transverse and the other longitudinal are excited in the spectral range ^ 410 — 440nm 
and ^ 600nm (only for e = 0.50) respectively. The longitudinal dipolar resonance is observed to blue shift 
with increasing aspect ratio. As the aspect ratio approaches £ ~ 1, the two transverse and longitudinal dipolar 
resonances merge. This can be attributed to the equal contribution of the respective polarizabilities. The relative 
intensities of the plasmon modes are also controlled by the strength of relative polarizability contributions. We 
must note here that the behavior of the resonances for nano-rods is a departure from that for ellipsoids where 
the transverse and longitudinal resonances get blue and red shifter respectively with increasing aspect ratio. ^ 
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Figure 1. The spectral scattering cross section Cscat for varying aspect ratios e exhibiting the quadrupolar and dipolar 
resonances. 

As mentioned earlier, the polarization characteristics described the polarimetric parameters like depolarization 
(A), diattenuation (d) and retardance (6) not only provide a wealth of information about the system under 
consideration, but like in plasmonic resonances, can be used as additional contrast parameters for high precision 
sensing applications. In the following, we investigate the spectral response of the polarimetric parameters in the 
context of their dependence on the aspect ratio. The calculations have been undertaken by preferentially orienting 
the long axis of the particle with the propagation direction of the incident radiation. The Mueller matrices 
obtained using the T-matrix method for preferentially oriented particles were subjected to the decomposition 
method described in the earlier section. 



0.8r 1.4r 0.8 




Wavelength X (nm) Wavelength X (nm) X (nm) 

(a) d at 6 = 50° (b) 5l at 9 = 50° (c) A at random orientation 



Figure 2. (a)Lincar retardance Sl and (b) diattenuation d with different aspect ratios for the nario-rod oriented at a zenith 
angle 0° to the laboratory frame of reference, (c) shows the net depolarization for orientation averaged, randomly oriented 
nano-rods. The scattering angle 6 in each case was kept to be 50° 

Linear diattenuation and retardance as mentioned earlier are the differential attenuation and phase difference 
between linear orthogonal polarizations. The linear diattenuation and retardance for silver nano-rods of EVSR 
20nm whose long axis is oriented along the propagation direction of the incident field, are shown in Figures 2(a) 
and (b). In 2(a), the linear diattenuation can be seen to peak at the plasmonic resonances while at the spectral 
overlap regions of two competing transverse and longitudinal modes is observed to show a sharp concavity or dip. 
However, it can be seen that as the aspect ratio e increases, the dip becomes broader. This can be attributed to 
the equal contributions of the transverse and longitudinal modes when the rod approaches a disk. In contrast 
to the diattenuation, the linear retardance spectrum can be seen to peak in regions where the plasmon modes 
overlap. Nano-rods as has been mentioned earlier have two polarization axes effecting the polarizability to have 
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two orthogonal components. The preferential orientation of the particle at /3 = 0° aligns the two polarization axes 
to coincide with the orthogonal polarizations of the incident electromagnetic radiation. This causes a differential 
attenuation of the two orthogonal polarization states while interacting with the metal nano-particle, the effect 
of which can be seen as higher diattenuation at the resonance bands. At the spectral overlap regions of the two 
bands, the polariz ability can be understood to be the vector sum of two orthogonal polarizations which would be 
commensurate with the low diattenuation as shown in Figure 2(a). However, at these spectral overlap regions, 
the phase difference between the two orthogonal polarizabilities is high which shows up as the increased linear 
retardance. 

Depolarization as mentioned before arises from multiple random scattering in the sample and as mentioned 
earlier is represented as the loss of polarization for 100% polarized incident electromagnetic radiation. In Figure 
2(c), the depolarization spcctriim is shown. The fact that the spectral signature of the depolarization closely 
resemble that of the linear retardance spectrum can be understood as the effect of the agglomeration of multiple 
retarders in random orientations giving rise to depolarizations. In fact, when the particle is preferentially oriented, 
the depolarization matrix Ma = Mjj. 

The diattenuation described in figure 2(a) was shown to exhibit a very sharp concavity at the spectral overlap 
region of the two competing transverse and longitudinal resonances. This sharp concavity as has been mentioned 
before is more prominent for aspect ratios much lower than 1 and the the; c;oncavity becomes spectrally wider as 
the aspect ratio approaches 1. This dip as we shall show presently, can be possibly used as a bio-sensor that can 
detect small fluctuations in the medium refractive index. Figure 3(a) represents the spectral sensitivity of the 
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Figure 3. Here in (a) Sensitivity of the diattenuation-dip of a nano-rod (e = 0.65, oriented at an azimuthal angle 30° 
with the laboratory frame) as a function of wavelength (A) is shown for different values of surrounding medium refractive 
index and (b) shows the figure of merit 7 a representative case. 

diattenuation dip of a silver nano-rod (e = 0.65, oriented at an azimuthal angle 30° with the laboratory frame) 
with refractive index of the medium. It can clearly be seen that as the medium refractive index rim increases, 
the diattenuation dip red shifts by ~ 3 — 4nm with every 10~^ increase in the n^. This spectral sensitivity 
can be understood while seen in conjunction with the scattering cross-section where, the spectral position of the 
longitudinal dipolar resonance also red shifts with increase in Um- The figure of merit (7) for the variation in 
spectral position of the diattenuation dip with varying medium refractive indices has been calculated as 



7: 
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and has been shown in figure 3(b). The figure of merit for the representative case shown in figure 3(a) is found 
to be 339.5809nm which signifies a spectral shift of 339.5809nTO for the diattenuation dip for unit change in the 
medium refractive index. 
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The fact that the (hattcnuation for plasmonic nano-particles is significantly higher than that of dielectric 
nano-particles of similar dimensions, and the spectral dependence of such diattenuation for metal nano-particles 
(arising from the spectral dependence of metal polarizabilities) can be exploited as a contrast mechanism for 
the detection of local refractive index variations in various inhomogeneous media where, the homogeneous vol- 
ume is significantly higher than the particle's volume. This diattenuation dip could consequently be used in 
various biomedical applications like real-time sensing of a biological cell's growth, determination of the sub- 
micron refractive index structure of tissues and provide various tools which could possibly used to probe disease 
progressions. 

4. CONCLUSION 

To summarize, characterization of LSPR in plasmon resonant silver nano-rods was carried out in this work using 
polar decomposition of scattering Mueller matrices in different configurations. The presence of quadrupolar 
resonance along with the dipolar resonance in the electromagnetic scattering spectra of such nano-rods show 
interesting effects on the spectral features of the polaximetric parameters. The decomposition analysis from 
the preferentially oriented particles revealed the presence of strong diattenuation (d) and linear retardance {S) 
effects in the regions corresponding to the resonances and in the overlapping regions of the resonances respectively. 
Distinct spectral features for both diattenuation and retardance were observed and are fomicl to be sensitive to 
the size and orientation of the nano-rods. These polarimetric parameters were further investigated for the nano- 
rods suspended in various dielectric media having different refractive indices and it was observed that small 
change in the mediuni refractive index leads to dramatic changes in the spectral response of both the retardance 
and diattenuation. However, the dip of the diattenuation arising between transverse and longitudinal dipolar 
modes was found to be maximally sensitive to the medium refractive index variations, with the figure of merit 
being 339.5809nm. This high sensitivity of the diattenuation dip, which experiences a spectral shift of more 
than SOOnm with a unit change in the medium refractive index has been further probed with an emphasis on 
biomedical applications. We have used the values of the medium refractive index corresponding to the local 
refractive index variations in biological systems and observed that the diattenuation dip does show sensitivity 
to such small changes and hence can successfully be exploited and used as the contrast mechanism in various 
biological applications. It is also worth noting that the unprecedented control over the polarimetric parameters 
for these silver nano-rods could possibly also be gainfully explotied for addressing fundamental questions in 
light-matter interaction, the details of which are under current investigation. 
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